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Because of the equal strength of the pairing potential mediated by one-gluon exchange for all 
partial waves to the leading order QCD running coupling constant and the nonlinearity of the 
gap equation, the non-spherical pairing in single flavor color superconductivity (CSC) can not be 
restricted in a single non-s-wave channel and the mixing among different angular momenta will 
occur. In this paper, we examine the angular momentum mixing in single flavor CSC with transverse 
pairing, in which the pairing quarks have opposite helicity. We flnd that the free energy of all non- 
spherical pairing states are lowered by angular momentum mixing compared with that contain 
p-wave only. But the amount of the free energy drop is numerically small. Consequently the most 
stable pairing state that respect the time reversal invariance remains the spherical CSL. 
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I. INTRODUCTION 



Color superconductivity(CSC) is expected to be the ground state in a quark matter of sufRciently high baryon 
density but low temperature P, [3, H, 0, H, 0, 0] ■ The most favorable pairing channel in the chiral limit is between 
quarks of different flavors. At the baryon density accessible in a compact star, however, the nonzero strange quark 
mass induces a significant mismatch between Fermi momenta of different flavors that suppresses the phase s pac e 
available for pairing. A lot of exotic CSC states or mixed states have been proposed in the literature [1, [§, [l^, [Tll.ll2j . 
But the result of their competition remains under debate. 

Since the single flavor pairing is free from the Fermi momentum mismatch, it has been considered as an interesting 
alternative in this circumstance. The corresponding gap function in this case should contain higher total angular 
momentum as required by Pauli principle, the obvious alternative is to consider the p-wave pairing(J = 1) jisl. fl4. fisl . 
m, [13) [ill ■ ^ large variety of pairing structures of the single flavor CSC has been discussed in the literature|l3l. Il4|. 
Except for the color-spin-locked(CSL) state, none of the energy gaps is spherical. All these works restrict the angular 
dependence of a non-spherical gap within a single partial wave, p-wave, which is valid if the pairing potential consists of 
the p-channel only, such as that in the superfluid ^He. The forward singularity of diquark scattering in QCD, however, 
brings in all angular momentum channels to the pairing potential. Because of the nonlinearity of the gap equation, all 
partial waves contribute to a non-spherical energy gap at the same order of the QCD running coupling. In a previous 
paper [l^. we examined this angular momentum mixing effect in the non-spherical CSC with longitudinal pairing, in 
which the two pairing quarks carry the same helicity. We found that the energy gap ( gap function evaluated at zero 
Matsubara energy ) takes the form 

where /i is the chemical potential, p is the direction of the relative momentum in a Cooper pair and the function /(p) 
is an inflnite series in spherical harmonics of all J's. In particular, we flnd 

/(p) = ei/3e"^[0.9866Pi(cos6l) - 0.0673P3(cos 61) + 0.0214P5(cos6') + ...] (2) 

in contrast with /(p) = e^/'^e~^Pi(cos 9) reported in the literature for the polar phase. The contributions from higher 
partial waves, however, are numerically small. The same effect is expected for the transverse pairing and the parallel 
analysis will be carried out in this paper. 
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The free energy density of the superconductivity with angular momentum mixing will be lowered relative to that 
containing p-wave pairing only. In the toy model considered in [19], which maximizes the forward singularity, the 
magnitude of the condensation energy is increased by a factor of 1.54. Without including the angular momentum 
mixing, the gaining factor of the magnitude of the condensation energy in CSC by forming the spherical CSL state 
compared with the non-spherical polar state is 1.54 for the longitudinal pairing and 1.14 for the transverse pairing [l^ . 
This raises the issue to recalculate the condensation energy of various non-spherical pairing states of the single-flavor 
CSC in the light of the angular momentum mixing and to compare the result with that of the spherical CSL. So we 
did in previous works [l9j [20l | for the polar and A phases of the longitudinal pairing. The formulation we developed 
can be carried over readily to the polar and A phases of the transverse pairing and the parallel analysis is made in 
this work. We found that, like the longitudinal pairing case, the amount of mixing is numerically small so the slightly 
reduced free energy of the polar and A phases is not low enough to compete with the spherical CSL phase. 

This paper is organized as follows. In Sec. II, we examine the quark propagator with transverse pairing. The 
condensation energy density and the gap equation will be obtained from the CJT effective action in Section III. The 
gap equation will be reduced to an nonlinear integral equation for the angular dependence with zero and nonzero 
azimuthal quantum number in Section IV and their numerical solutions will be given there as well. We conclude in 
Sec. V. Some technical details about the full quark propagators and the relevant spherical harmonics are deferred to 
Appendices A and B. Our convention for the metric tensor is g^^^ = diag{l, —1, —1, —1}. Our units are h = c = = i 
and 4- vectors are denoted by capital letters, K = = (w, k) with i/ the Matsubara energy, which becomes continuous 
at zero temperature, T — 0. 



II. THE QUARK PROPAGATOR FOR TRANSVERSE PAIRING 

Though the general structure of the quark propagator in CSC state has been discussed extensively in the literature 
[2l[, the formulation we developed below is particularly suitable to explore the angular momentum mixing in a 
transverse pairing state. 

We start with the thermal Nambu-Gorkov quark propagator S{X), defined as 

where and ^ are the Nambu-Gorkov(NG) spinors, P — [iv, p), P X = — i^r — p-x with v the Matsubara energy and 
T the Euclidean time, T imposes the Euclidean time ordering and (a, (3) run over the NG, Dirac and color indexes. 
At zero temperature, = The NG spinors are defined by 

* = (V;,^c), (4) 

where, "0 and tp are written in terms of the ordinary quark spinors, ipc = Cip and ipc — ipC arc the charge-conjugate 
spinors with the tilde standing for the transpose, and C = 17^7°. It follows from the definition ^ and the self- 
conjugation relation ^! = C^** with 

C - ( 1 ) ■ (5) 



172 
that 

S*{P) =C^S{P)C (6) 

with P = {iv, — p) and 

S{P) = ToCS{-P)roC (7) 

with 

The relation © was proved in 0| and the proof of ([7]) proceeds similarly. The NG structure of the inverse propagator 
is 



, f K+{P) $-(P) \ 



(9) 
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where each block is a 4 x 4 matrix with respect to Dirac indexes and a 3 x 3 matrix in color indexes. We shall focus 
on the ofF-diagonal blocks <i>*(P) below. 

It is convenient to work in the chiral representation where 



V'= ( ), ^ = (10) 



and their charge-conjugate counterparts 



^c=( 'fl}.), ^Pc = {iXR<^^-tXL'y'), (11) 

with xl{xr) the two-component left-handed(right-handed spinors) and the second Pauli matrix. The matrix 
(P) are decomposed accordingly. For example, 

with each block a 2 x 2 matrix with respect to spinor indexes. For transverse pairing between different chiralities, 
$^_i_(P) — <I>^_(P) = 0. Next, we consider the off-diagonal block (^^_(P) and write 

$;_(P) = ip{P) + <^o{P)a- ■ p + ^+{P)a- ■ e+(p) + ip-iP)o- ■ e_(p), (13) 

where the circular basis e±(p) — -^{Oztitp) with 6 and (p the unit vectors of the spherical coordinates of p. We have 

p = (sin cos iy9, sin0sin(/3, cos ^?) 
— {cos 9 cos (f, cos 9 simp, — sin 9) 

(fi — (— sin(p,cos(p,0), (14) 



e+(-p)-e_(p) e+(p)* =e_(p), (15) 

and 

P X e±(p) = Tje±(p). (16) 

Performing a rotation about p by an angle 7, wc find 

$;_(F) ^ e-^^'^-P$;_(P)e*^'^-P ^ p{P) + MP)o- ■ P + e-^T^+(P)(T • e+(p) + eJ''ip^{P)(T ■ e_(p). (17) 

For a Copper pair of two quarks or two antiquarks of opposite helicity, the projection of its angular momentum in 
the direction p has to be ±1, which implies that (p{P) — Pq{P) — 0. Up to now, the coefficients (p^s remain 3x3 
antisymmetric matrices in color space. For the polar phase and A phase considered in this work, we set 

$;_(P) = o--*(P)A2 (18) 

where 



*(P) = MP)e+{p) + </.-(P)e_(p). (19) 

and A2 is the second Gell-Mann matrix. 

The Dirac structure of other off-diagonal elements in Eq.® can be deduced from the relations ^ and ([7]). We 
now could write down the inverse full quark propagator in chiral representation using the Dyson-Schwinger equation 
and ignoring the contribution from the wave-function renormalization, and reads 



iiy + jj. — (T ■ p 

i + (T • p — A2(T • 

_ A2cr-**(P) 

-A2cr-**(-P) iiy-^ + cr-p 



o-VpN W + ^l + CT-p -A2(T-*(-P) I 

^ ^^)-\ A2cr.**(P) ,-w_„_^.r. I l^tJi 
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Sii 


S12 


Sl3 


5*14 


S21 


S22 


S23 


5*24 


S31 


S32 


S33 


5*34 


S41 


S42 


S43 


5*44 



Inverting the matrix ([20]) (see the appendix A for details) , we find that 

/ Sii S12 Si3 5*14 

QfP\ — I "^21 S22 S23 5*24 I ^rjiN, 
'^K^) - I Q„ Q„„ Q„„ Q„. I l^-LJ 

with 

Sii = Si3 = S22 = S24 = S31 = S33 — S42 = S44 = (22a) 

a ^ iiy-{p + fi) l + fT-p ii^+{p- fi) 1 crp 

u^ + {p + ^if + 2\l^^{~PW_{~P) 2 ^^ + {p-iiY + 2\l<j,+ {-P)cpX{-P) 2 ^ ' 

o ^ AaCT ■ e+(p)(?!)_(-P) ^ 1 - <t ■ p Aacr ■ e_(p)0+(-P) ^ 1 + cr ■ p 

" ly^ + ip + ^i)^ + 2Xl^-i~P)r-{-P) 2 z/2 + (p_^,)2 + 2A2</)+(-P)0;(-P) 2 ^ 

a ^ z;/ - (p + fi) lap iu+{p- ^) 1 + q- p 

I/2 + (p+^)2+2A20_(P),/.^(P) 2 i.2 + (p_^)2+2A20+(P)0;(P) 2 ^ ^ 

^ A2T ■e-(p)0-(P) ^ l + (T-p A2cr ■e+(p)(?!)+(P) ^ 1 - cr ■ p 

i/2 + (p + ^)2+2A20_(P)01(P) 2 ^ z.2 + (p_^)2+2A20+(P)0;(P) 2 ^ ^ 

^ ^ A2cr-e_(p)C(^^) l + cr-p A2cr ■ e+(p)0;(-P) 1 - cr ■ p 

i.2 + (p + ^)2 + 2A20_(-P)01(~P) 2 z/2 + (p-/x)2 + 2A20+(-P)0;(~P) 2 ^ ^ 

c ^ iiy+{p + fi) lap iu-{p- fi) l + cr p 

l/2 + (p + ^)2+2A20_(-P)01(-P) 2 z/2 + (p_^,)2+2A20+(-P)0;(-P) 2 ^ 

^ ^ A2<T.e+(p)01(P) l-<T.p A2<T.e_(p)0;(P) 1 + ,t ■ p 

l/2 + (p + ^)2+2A20_(P)01(P) 2 ^ i/2 + (p_^)2+2A20+(P)0;(P) 2 ^ ^' 

c ^ z;/ + (p + /i) 1 + ■ p it/ - (p - /i) 1 - o- ■ p , . 

i/2 + (p + ^)2 + 2A2</,_(P),/.l(P) 2 t.2 + (p_^)2 + 2A20+ (P)0;(P) 2 ■ ^ 

We notice that associates with particles and (j)- with antiparticles. The is because we designate implicitly the 
momentum p to the quark of positive helicity in a Cooper pair. Since the contribution from the Cooper pairs of 
antiquarks goes beyond the subleading order, we shall approximate (f>-{P) ~ and introduce the gap function (j)(v, p) 
via 

MP)^^^i^,P)- (23) 

As is the same with 2SC, the solution to the gap equation gives rise to an even function in the Matsubara energy, we 
have M-P) = -^<^{'^,-p), MP) = -^M^P) and <^+(-P) - -^M,-p). 



III. THE GAP EQUATION 



In this section, we shall derive the gap equation for the transverse pairing. The starting point is the CJT free 
energy, which can be approximated in weak coupling as[i, 

r = r„ + QF (24) 

where, r„ is the free energy of the normal phase and F the condensation energy density, which is the part of T 
responsible to the gap equation. The CJT formula for F reads 

F = -^[TrlnP>-i - TrlnD,;! - TrlnS'-^ + Trln5(7^ - Tr(5(7^S' - 1)] (25) 

where D and S are the full propagators for gluon and quark respectively. So is the tree level quark propagator and 
Dn is the hard-dense-loop (HDL) gluon propagator. Following the procedure of 0, [111, we approximate 

TrlnD^i - Trln[i:>^i] ~ TT[DnSU] (26) 
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where 



(5n = n - n„ 



(27) 



with n„ the HDL rcsumnicd ghion self-energy in normal phase. Eg. ((26)) represents the Hartree-Fock energy of the 
mean-field theory. In the Coulomb gauge, the HDL gluon propagator is 



A.,oo(^) - Di{K), Dr,,o^{K) = Amo(^) = 0, = ^ kk)Dt{K) 



(28) 



where Di^t are the longitudinal and transverse propagators respectively and are diagonal in adjoint color space, 
i.e. \ — S°'^Di,t. Consequently, we only need the 00-component, H""(i4r), and the transverse projection of the 
ij-components. 



The gluon self-energy in super phase reads 

p 



(29) 



(30) 



where P' = P - X, ^ T / (C^ and 



/ a^Ta \ 

a^'Ta 

-cr^fa 

V -<J^^fa J 



(31) 



with a^' = (1,<t), ct'' = (1, -cr), and Tq = Aa/2 (a = 1, ...,8) is the ath SU{3) generator with Xa the ath Gell-Mann 
matrix. Substituting Eq. ipTjl and pTjl into Eg . ([501) . we have 

KbiK) = ^9'Y.^T[a^TaS2i{P)a''nS2iiP')+a^TaS,2iP)a''nSMP') 
p 

-a^TaS23{P)a''nSii{P') - a^TaSu{P)a''nS^2{P') 
-a^faS^iiPyTAsiP') ~ a^fa532(P)a"Tb5i4(P') 



+<J^'faS43{P)<j''ftS43{P') + a^^faS3i{P)&''nS3i{P') 

Substituting Eq.jlll) to Eq.(l32l), we find that 



(32) 



Tr 
Tr 
Tr 
Tr 
Tr 
Tr 
Tr 
Tr 



<T^'TaS2l{P)a^nS2l{P') 

ij''TaSi2{P)a''nSi2{P') 

a''TaS23{P)<j''nS4l{P') 
<J^TaSii{P)<j''nS32{P') 
<j'^faS4l{P)<j''nS23{P') 

a^faS32{P)^''nSu{P') 

<J>'faS43{P)<j''nS43{P') 
nS34{P)a''nS34{P') 



= Tr 



Tr 
1 



iTr 
2 

iTr 
2 



wi(p,p'), 
^(P^P'); 

a^<T . e+(p)0(i.,p)i^|-P(7V . e_(p')0*(i^', P')^^|-^ 
a'V . e_(p)</.(z., -p)l±^a-a ■ p+{p'W{i^\ -P')^— ^ 



a^<T . e_(p)(/.*(j.,p)i±^a''^ • e+ip'^i.', p')^-^ 



a^<T ■ e+(p)0*(^., -p)L^-^a'^a ■ e_(pW, ^p') ^ + ^ P' 



Tr 
Tr 



ul-cr-P i.l-cr-p' 
c ;: cr 

2 2 

_ul + cr-p_^l + cr-p' 



2 

W'4(P,P'): 
W4(P,P') 



W2(P,P')) 

U'2(P.P'): 
W3(P,P'). 



(33) 
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The trace over color space can be readily performed and we have 

( 1 iiy±£, ii^'±£,' 



wia{p,p') = w'l 4(-p, -p') = < 



2 Z/2 + ^2 ^ A2 I/'2 + £'2 + A'2 ' 

1 r 1 w'±i 1 



il. ip £ jy'2 + £'2 _^ A'2 =F + ^2 + A2 



1 iiy±£, 



iv' ± 



1 1 



I 6 ^2 ^ ^2 A2 ;y'2 + f 2 _^ 3 =p £ ii,' =p £' ' 



a = 6 = 1,2,3 
a = 4, ...,7 

a = = 8 



(34) 



and 



W2(P, P') = W2(-P. -P') = W^3(P, P') = W3(-p, -p') = "( 0' 



2 1/2 + £2 ^ ^2 j,/2 _^ £.2 ^ /^.2 ^ 



1 



1 



6 1/2 + + A2 j/'2 + 



a = 6 = 1,2,3 

a = 6 = 4,.. .,7 (35) 

a = 6 = 8 



where we have defined ^ = p — n, A'^ = |0(i/, p)P- The ± on the RHS of Eq. ([M|) correspond to wi and W4 respectively. 
In weak coupling, the dominant contribution arise from the quasi-particles, therefore we have ignored that from the 
quasi-antiparticles. The trace over Pauli matrices on the RHS of Eg. ([55)1 are evaluated. We have 



ol±(T-p ol±cr-p' 

Tr^i^^a^i^^] 

L o O J 



Tr[."i^.°i^]4(l+p.p') 
= Tr[.^i±fP.^i±fP:]4(3-p.p') 



kikjTila 



l±(T-p ,l±cr-p' 



It follows from Eq.((Tl]) that 



kikjTila' 



l±0-p_,l±(T p' 



= i[l-p.p' + 2k.pk.p']. 



• e±(p) = e±(p). 



and the trace in Eq. (|33p is simplified. The terms we need are 

, ,,l+(T.p' 



Tr[aV.e+(p)0(i/,p) 



l^-PaV.e_(p')0*(^',p')- 



Tr[c7 V . e+ (p)</)(i/, p) V • e_ (p')0* (i/', p') ^^T^' ] 



20(;.,p)0*(z/',p')e+(p).e_(p') 
^2</)(j/,p)0*(^',p')e+(p).e_(p') 



(36) 
(37) 
(38) 

(39) 

(40) 
(41) 



and 



^,%TrKT • e+(p)0(t/,p) ^ r^ ^^^ • e-(p')</>*(^',p')^^^] 



4(pp' - p • p') 



(P P 



,M2 



z/,p)0*(z.',p')e+(p) • e_(p') - 20(i/,p)r (^^',p')e+(p) • e_(p') 



(42) 



The other cases in Eq. (|33p can be evaluated along the same lines and we do not write them down explicitly here. 
We now turn to the Hartree-Fock energy density in Eq. (|26p and divide it into two parts 



HF 



1 

m 



where the first term comes from magnetic gluons and the second term from electric gluons and we have 

= ^19' E A(if)(% - kk,)TT[Tis{p)nsiP') 

P,P' 

= -^ff' E Di{K)Ti-[rlS{P)rlS{P') 

P,P' 



(43) 

(44a) 
(44b) 



As we discussed in our previous paper [19j, the net contribution from Eq. ([34)) and that from the diagonal block 
of the quark self-energy matrix is to add a wave- function renormalization factor in the gap equation(26l. [27l . [28l . [29| . 
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which is subleading order. It will not interfere with the angular dependence of the gap function as we have already 
seen in our previous paper [T9|. Therefore, the contribution from Eq. (|34p will be ignored and only that from Ea.(|35p 
will be considered here. We find 

r ^ _ ^ 2 V- r) / y^ PP'-P- P' '/'(^: P)(t>* {v'; POe+ (p) ■ e_ (pQ + c.c. 

3^ *^ ' (p-p')2 +^2 + A2)(iy'2 +^2 + A/2) 

+ {a term with p — p, p' — p'} (45) 



and 



= E ^^(^) ( J + ^2 + A2K.-2 V + A-2) + term with p -> -p, p' - -p'} (46) 



In weak coupling, the approximation p ~ p' ~ is justified, and we can perform the integral over ^ by picking up 
the the poles in the diquark propagator. Assuming further that the gap function depends only on the energy and the 
orientation of the momentum, we find 



2 4 

HF = -77.9 M 



1 2 i f dv f dv' f (Pp f (Pp 



6-^ 7 27r 7 27T J (27r)2 J {2tt) 

In order to evaluate the 3rd and 4th terms on the RHS of Eq.(l25|). we employ the identities TrlnS*^^ = lnDct5^^ 
and 

Det ^ ^ ) = Dct(-CB + CAC-^D) = Det(-BC + BDB-^A) (48) 

with A, B, C and D arbitrary matrices and B,C invertible. We have 

1 /m 0-1 rr^ 0-1 \ n 2 f dv I' (Pp 



^^-(Trln5- - Trln5-) = J - V^^T^) (49) 

The evaluation of the last term on the RHS of Eq. ([^S)) can be readily performed by using the full quark propagator 
and we end up with 

^2^P l^^l ^ (50) 

2ri ^ ° ' i 27r 7 (2^)2 Vj^2 + A2 ^ ' 

Therefore, the condensation energy density reads 

f f f ^2^ f M^f,,,.. ..f ^ ^„0(i^,p)0*(i^',p')e+(p) •e_(p') +C.C. 



- / dv / dv' / d'p / d'p'V{v-v',p-p 



647r4 J J J J ' '^i-/ _^ A2Vi^'2 + A'2 

+ iJ/../.^p(M-^) (51) 

where, the pairing potential 1^ contains the contributions from both magnetic and electric gluons and reads 

V{v -v',p- p') =Dt{v- v', p - p') + Di{v ^v\p- p') (52) 

with p — p' = ^i{p — p')- In the energy-momentum region for Cooper pairing, \v — v'\ << |p — p'|, Dt is dominated 
by Landau damping and Di by Debye screening. We have 

V{v^v',p-p') = . n- + 7 ^ r (53) 

(p-p')2+>i{^ ip-p'y+< 

where the Debye mass square mj-, = —^^pr—- 

To obtain the gap equation, we introduce the quantity, 

C = -^^= (54) 
^ Vv^ + A2 ^ ' 
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The minimization of the condensation energy density with respect to ^ gives rise to the gap equation 



3gV 

47r 



,,</>(z.',p')e+(p')-e^(p) 



+ A'- 



(55) 



The forward singularity corresponds to the divergence of the angular integration when the coupling g in the denomi- 
nators of is set to zero. To isolate out the singularity, we write the integral on RHS of ([55]) as 



du' J d'^p'V{v-v',p-p') 
dv' [ d^p'V{iy-iy',p-p') 



^iy'^ + \^iiy\pW 
(i.',p')e+(p')-e-(p) 



(56) 



The coupling g in the second integral can be set to zero without introducing divergence. So we did and the gap 
equation reduces to 



dv'{ -In 



2 ^v'^ + \ct>{v',pW 



+ 



47r 



•2-1 



d^p 



1 - p • p' 



w.p) 



0(^^^pOe+(pO■e_(p) ^ 

v^^'2 + i0(i.',p')p v^t^trt;^ 



(57) 



where we have carried out the p' integral of the first term of (|56p and have inintroduced an infrared cutoff ujq ~ g/i. 



IV. THE INTEGRAL EQUATION FOR ANGULAR DEPENDENCE AND ITS SOLUTION 



In the expansion Eq. (ll9|) . the circular polarization basis e± carry singularity either at the north pole or at the south 
pole of the spherical coordinates. The coefhcicnt functions that remove this singularity are the Wigner D-functions, 
instead of the ordinary spherical harmonics. The Winger D-function corresponds to the irreducible representation of 
the rotation group 30] and is defined by 

D'Lm' (a, /3, 7) = exp(-iTOa)rf;^^, (^)exp(-M7i'7) (58) 

with 

Jm\e~''y^\Jm' > . (59) 

Its complex conjugate, £*„*„/(«, /9, 7), is the wave function of a symmetric top of definite J^, and J • with C, the 
symmetry axis. A diquark pair of opposite helicity of a definite angular momentum represents a parallel case with 
C of the symmetric top replaced by direction of the relative momentum of the pair, p. Since J • p = ±1, the wave 
function of a transverse diquark pair ought to be proportional to Z3^*j_]^. Here, however, there is a question when 
carrying this wave function to the diquark pair. A general rotation is defined by three Euler angles a,/?, 7 but the 
gap function depends on only p, which is defined by the two polar angles (f and 9. In what follows, we shall identify 
the first two Euler angles with the polar angles, i.e. a = ip, P = 9, and set 7 = 0. As is shown in the appendix B, 
the real and imaginary parts of D'J*^^{ip, 9, 0)e±(p) are proportional to the vector spherical harmonics, which are free 
from the singularities at 6* = and 9 = tt. Therefore the partial wave basis of the gap function (|23p is D'^i{(p, 9, 0). 
The partial wave amplitude of the pairing potential, V{i' — ly' ,p ■ p')e+(p') • e_(p), is given by the integral 

J d'p J d'p'Di^,{^,9,0)V{,y~,y\p-p')e+{p')-e^{p)Di*^,{^,9,0) 

= Jjfj J d'p J d'p'V{.-.',p-p')Y^D'U^,9,0)Di\i^,9,0) 

= ^ / d'p J dVn---',P-p')dfi(e)cos2| 

= 2JTT'^"''"'/ d cos QV{iy-iy', cos Q)dii{e) cos'' ^ (60) 
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with cosQ = p • p', which corresponds to expansion coefficients of V{iy — i^', cos 0) cos^ ^ in the series of d(i{Qys. 
The first equaUty of ([SO]) follows from the Wigner-Eckart theorem and the second one from the addition formula 
(B13) (proved in Appendix B) 

J2 D'ZM 0, 0\ 0)e+(p) • e_(p') = di^{Q) cos^ |. (61) 



We have[2a 1271,1231 

V{v - p - p') cos^ I = Y.i2J + l)di,{&) + ^ £(2 J + l)c'jdi,{Q) (62) 

where, lUc = "'"'^^'^^^ ^ and c'j is given by 

forJ^l 

dcose ^"7^^"" ^ ^ = <; ^1 1 (63) 



1 



1 - cos e 1 _ I t for / > 1 

I Z^n + 2J(J+1)' *°'^'^>^ 



Unlike the magnitude of the angular momentum, J, that is mixed in the solution because of the nonlinearity of the 
gap equation. The azimuthal quantum number, — m, remains conserved. For a gap function of a definite to, 

(/.(:/,p) = Vm('^,cos0)e™'^ (64) 

with ?/;„(;/, cos 0) — real, the gap equation (j57p implies the following integral equation satisfied by V'm('^, cos0) 



1 ( 1 

V'm(i^,a;)= -.9^y di^'|-(ln 



1 

dx'- 



1 \Im{x,x')'4}„i(v' ,x') i'{v',x) 



\x — 



(65) 



where a; = cos 0, a;' = cos Q' and Im{x^ x') =\x — x'\Im.{x^ x') with 



, . 1 e+(pO-e_(p)e-('^'-^) 

I„i{x,x)^— dip : — . (66) 

2tt J„ 1 - p • p' 

Substituting the explicit expression of e+(p') • e_(p), we find that Im{x,x') is real and 

I-rnix,x') = I,n{-X,-x'). (67) 

The reduction of the gap equation is parallel to that for the longitudinal pairing case. We shall outline only the 
main steps here and refer the interested reader to the details in the section IV of [19]. We start with the approximation 
in [13 

In^^ ^ (68) 

W-i^l w>\ 

with t/> = max(j/, v') and introducing the new variable 

e = ln^. (69) 

Twice differentiation of Eq. (|65p with respect to ^ yields an ordinary differential equation 

d'^tpm{^,x) g'^ip^{^,x) 



together with the boundary condition 



(70) 



dll^mj^x) ^ ^ 
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as ^ — > oo, valid for all x. The solution of Eq. ((70| up to the subleading order reads: 

V'm(C, x) = Ao/m(x) x)u{£,, x) + B{£,, x)v{^, x) 

where 



COS g[b(x) — for ^ < b{x) 
1, for £_ > b{x) 



v{^,x) 



— sin g[b{x) — ^] , for ^ < b{x) 
gi, for C > b{x) 



(72) 

(73) 
(74) 



and 



^l + u2(f,x)e2l«'-bWI 



(75a) 



v/l + u2(f,a;)e2l«'-K=^)l 



^'(^',x). 



Here, we have introduced the energy gap of 2SC 



2048\/27r4/i 



(75b) 



(76) 



and the quantity b{x) — In j ''^(■^■)| ■ The angle dependent function fm{x) will be determined below. We notice 
that at zero Matsubara energy, ^ ^ oo, ^pm{0,x) — Aofm{x). Substituting Eg. ([72)1 back into Eq. (|65|) . we derive the 
integral equation for fm{x) 



/™(x)ln|/m(x)| + 3 / dx- - 

_i \x — x'\ 



fm{x) X^(x,X ^fni{x ) 



0. 



(77) 



The solution of (I77p contains all partial waves with J > max(l, |m|) and we shall concentrate on most favored pairing 
channels, m = 0, ±1. In '^He, m = refers to the polar phase and m = ±1 to the A phase. We have 



Io{x,x') 



XX 



1 



■( 



VI - x'"^ vT^^ 



2Vl^^VT^^ 2\x-x'\\^T^ 



and 



I±i{x,x') = - 



1 



l±x' l±x 



X — x'\\ 1 ± X 1 ± x' 



1 



1 



l±x l±x' 



(78) 



(79) 



Because of the symmetry (1671) . given a solution of the integral equation of the azimuthal number m, fm{x), the 
function fm{~x) solves the equation of the azimuthal number — m. For the solution starting with the J = 1 harmonics, 
we expect that 



f-m{x) = fmi-x) 



(80) 



SO fo{x) will be an even function of x. Since dQi{d) is proportional to the associate Legendre function P]{x), only 
the partial waves of odd J contribute to /o(a;) and the higher harmonics starts from J = 3. As to f±i{x), all higher 
partial waves with J > 2 will be mixed in, so the free energy reduction is expected to be larger than that of fo{x). 

As we have done in [l^, this type of integral equation can be solved by a variational method. Upon the leading 
order solution of the gap function and the definition of angular dependent factor, the condensation energy density 
Eq. ([5T|) turns out to be a functional of /(x)(see appendix B in [19]). we have 



F : 



(81) 
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with 



-F[/„.(x)]= l'^dxfU^)[\n\f„,{x)\-^] + ^ j'^dx l'^dx'j-^^[f^{x)-2U{x)U{x')I^{x,x') + fUx')} (82) 

It can be easily verified that the variational minimization of Eq. (|5^ does solve Eg. ([77)1 . Before the numerical calcu- 
lations, we consider a trial function for each case, to = 0, ±1. For m — 0, the trial function is 



fo{x)=codUe)^co\i^-^ (83) 



and substitute it into the condensation energy density Eg. ([5^ . The minimization yields 

1 

V2 



Co = -^eee 2 (84) 



at which 

-1.214 X 10""* (85) 

For TO = ±1, the trial function is 

/±i(x)=c±4ii(cos0)=c±ii^ (86) 

The same procedure will give rise to 

c±=e*e"5 (87) 

at which 

-8.084 X 10"^ 



The target functional corresponding to to = ±1 has the same value, which is also true when the situation goes to 
angular momentum mixing. The condensations with to = ±1, therefore, have the equal priority when the temperature 
cools down to the critical temperature. 

The trial function E g. (1 831) and Eg. (j86p correspond to what people carried over from the polar phase and A phase 
of '^He respectively [l^ Il4|. which contain p-wave only. The corresponding condensation energy are not optimal as we 
shall see in the following, which will be lowered further by including higher partial waves. 

Following the same procedure in section five in 19!|, we can obtain the numerical solutions of Eg. ([77)1 and the 
corresponding condensation energy density. In Fig[l] we show the angular dependence of the gap function for the case 
to = 0, the dashed line and the solid one are the trial function and the numerical solution to Eg. (|77p respectively. 
They are depart from each other slightly, however, the condensation energy density corresponding to each case has 
qualitative difference. We find in this case the minimum of target functional 

T = -1.236 X 10"" (89) 

which drop from Eg. ([55)1 by 1.81 percent. 

The same is true in the case to = ±1. In this case, the trial fimctions (dashed lines) and the numerical solutions to 
Eq. ([77)1 (solid lines) have been plotted in FigHj We find the minimum of the target functional corresponding to solid 
lines 

T = -8.494 X 10"^ (90) 

which drop from Eg. ([55)1 by 5.07 percent. 

As we have found in [19|, the drop of the condensation energy with longitudinal pairing by angular momentum 
mixing is a small amount. Here, with transverse pairing, the situation is also true. With these small amount drop in 
condensation energy, the non-spherical pairing of the polar and A-phases in single flavor CSC remains not energetically 
more favored than the spherical pairing state CSL. 

It is instructive to examine the angular momentum content of our solutions according to 

oo 

/™(cos0) = J2i2J+l)bjdiM (91) 
J=l 
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FIG. 1: The angular dependence of the gap function for m=0. The dashed and the sohd hues correspond to trial function 
Ea. (|83|l and the numerical solution to Eq.((77} respectively. 



0.016 
0.014 
0.012 
0.010 
0.008 
0.006 
0.004 
0.002 
0.000 
-0.002 



m=1 



m=-1 



I ^ I ^ I ^ I ^ L 

-1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 



FIG. 2: The angular dependence of the gap function for m=±l. The dashed and the solid lines correspond to trial function 
Ea. (|86p and the numerical solution to Ea. (|77|l respectively. 



for m = 0, ±1, wc found 



and 



/o(cos 61) = ^e5/^e~^/2^1.00854i(6i) + O.O426<(0) 
v2 



(92) 



/±i(cos6') = ei/='e-^/2[1.0269diii(6') T O.O9434ii(0) 



(93) 
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TABLE I: The condensation energy of various non-spherical pairing states in unit of the transverse CSL condensation energy 
with or without angular momentum mixing. 





Without mixing 




With mixing 




F/FcsL 


Longitudinal 


Transverse 


Longitudinal 


Transverse 


Polar 


0.65e-^7e-^ 


0.88 


0.67e-^Ve-« 


0.9 


A 


0.88e-iVe"' 


0.65 


0.9e-iVe-' 


0.68 



V. SUMMARY AND DISCUSSIONS 



If the pairing potential extends to more than one partial waves, a non-spherical energy gap cannot be restricted 
within a single angular momentum channel because of the nonlinearity of the gap equation. The corresponding free 
energy will be brought down by the mixing. In case of the color superconductivity mediated by one-gluon exchange, 
the pairing potential mediated by one-gluon exchange for all partial waves equal in leading order in QCD running 
coupling constant. The angular momentum mixing will occur to the leading order of the gap function and free energy 
of all non-spherical pairing states in the literature needs to be recalculated. So we examined the angular momentum 
mixing in this paper for transverse pairing and in a previous paper for longitudinal pairing. We derived the nonlinear 
integral equation for the angular dependence of the gap function and solved it numerically. We found that the effect 
is numerically small, which may be attributed to the descending pairing strength of higher partial waves beyond the 
leading order. Our results on the condensate energy of various non-spherical pairing states of single quark flavor are 
listed in the table I below, in which the ratios without mixing are also listed for comparison. 

While the free energy of all non-spherical pairing states are lowered by the angular momentum mixing, the amount 
of the drop is too small to beat the spherical CSL state of transverse pairing. The energy balance between the CSL 
and non-spherical states may be offset by anisotropy or the violation of time reversal invariance. In a compact star, 
the magnetic field and the stellar revolution provides such a circumstance. The magnetic field, however, is unlikely 
to discriminate the non-spherical pairing states from the spherical ones in the bulk because of the Meissner effect in 
the CSC phase [l^. In addition, it suppresses all single flavor pairing because of the critical field. As to the stellar 
revolution, the optimal partition of the angular momentum into the motion of the center of mass and the spin of 
Cooper pairs remains to be determined. But the non-spherical condensate will response differently from the spherical 
ones because of the preferred direction in the former. 

The slightly complicated non-spherical CSC phase, the transverse planar phase exhibits high potential to be able 
to compete CSL phase since the gain of the condensation energy from the former to the latter is only 2 percent . 
falling within the range of the percentage increment of the condensation energy because of the angular momentum 
mixing discovered so far. Technically, the color structure of the order parameter involves two antisymmetric Gell- 
Mann matrices and the present formulation needs to be generalized. This is currently investigated and the result will 
be reported shortly. 
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APPENDIX A: THE ENTRIES OF THE FULL QUARK PROPAGATOR 



In order to obtain the full quark propagator, we introduce 



70 




ro= ( [;\J and V= I ^ , ^ ^ 1 (Al) 
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Using these two matrices, we can transform the inverse quark propagator into a block diagonal form 

G{p) = vs-\p)roV-^ = 

where the matrix T is given by 



T-\P,X2) 

T-\-P,-X2) 



T{P,X2) = 



iiy + fi-cr-^p A2(T • 
X2(T ■ ^* (P) iv ~ 11 — cr ■ ■p 



(A2) 



(A3) 



and can be inverted readily by using the projective operators P± — (1 ± cr • p)/2. Then we have the full quark 
propagator 



S{P) = ToV-^G-^{P)V 



I S\\ S\2 5*13 5i4 \ 

5*21 ^22 '5'23 '524 1 

5*31 ^32 5*33 5'34 I 

\ 5*41 542 'S'43 'S'44 / 



with each entry given by 



5ll — 5l3 — 522 — 524 — 531 — 533 — 542 — 544 — 

iv — e{p + efi) 1 + ea ■ p 



e± 



{p + e^l)^ + Xla- ■^{-P)(T ■^*{-P) 2 



5i4 = -A2<T.VI/(-P)^. 



1 



1 + ecr • p 



e± 



(p-e^)2 + A2cr-**(-P)cr-*(-P) 2 



e± 

S23 = X2a-<f{P)J2 
532 = -A2<T.**(-P)^ 



iiy + e{p — e/i) 



1 + ecr • p 



i/2 + (p-e/^)2 + A2cr-*(P)(T-**(P) 2 

1 1 + ecr • p 

+ (p + e/i)2 + Xja- ■ ** (P)cr • *(P) 2 

1 1 + ecr • p 



(p + e^)2 + A2cr-*(-P)cr-**(-P) 2 



534 = 



w + e(p — e/j,) 



1 + ecr • p 



e± 



'2 + (|5-e^)2 + A2cr-**(-P)cr-*(-P) 2 



541 = A2^T.**(P)^ 

e± 

543 = - E 



1 + ecr • p 



e± 1^2 (p _ e//)2 _^ x^a- ■ *(P)cr • (P) 2 



w + e(p + e^) 



1 + ecr • p 



e± 



i/2 + (p + e/x)2 + A2cr-**(P)cr-*(P) 2 



(A4) 

(A5a) 
(A5b) 

(A5c) 

(A5d) 

(A5e) 

(A5f) 

(A5g) 

(A5h) 

(A5i) 



In order to get rid of the Pauli matrices in the denominators of the propagator in Eq. (|A5p , we expanded ^ in terms 
of the two circular polarization basis as Eq. . We have 



a ■'>lf{P)(T ■^*{P) = *(P) • **(P) +icr • [^(P) X **(P)] 

= (1 + • p)0+(P)0;(P) + (1 - ^ • p)0_(P)C(J^) 

Therefore, the entry Eqs. (jA5dp can be written as 

iv - {p + Ij) 1 - cr • p 



(A6) 



521 — — 



u^ + {p + ^iY + A2[(l + ^ . p)0+(P)0;(P) + (1 - o- • p)(^_(P)r-(P)] 2 



iv+[p- /i) 



1 + cr p 



^.2 + (p - ^)2 + A2 [(1 + ^ . p)</,+ (P)0:; (P) + (1 - <T • p)(/._ (P)^l (P)] 2 



It can be shown that 



l±cr-p l±cr-p 
/('^ ■ P) ^ = /(±1) 7^ 



(A7) 



(A8) 



15 



with f{(T ■ p) is an arbitrary function of cr • p as long as /(±1) is well-defined. Using this identity, we find 

c ^ iiy~{p + n) l-a j, iv+{p- n) 1 + ct ■ p 

i/2 + (p + ^)2 + 2A20_(P)01(P) 2 + {p - fiy + 2XlMP)KiP) 2 ^ ' 

The other entries of the quark propagator can be simplified along the same lines. Furthermore the entries with ^ 
or ^* in the numerators can be reduced following the identity ([39]). The final expression of the quark propagator is 
presented at the end of Sec. II and is used to calculate the condensation energy density in Eq. (|^ . 



APPENDIX B: VECTOR SPHERICAL HARMONICS AND ADDITION FORMULA 

To each momentum p, one defines two transverse spherical harmonics (sij : 



Y«(P) = X,„,(p) = _— Ly,,„(p) (Bl) 
and 

Yf2(p)=pxX,„,(p) (B2) 

where the angular momentum operator L = —ip x Vp with Vp the gradient operator with respect to p. It follows 
from the relation bewteen the scalar spherical harmonics and the Wigner D-function that 



1 



^^"^^^ = V Jj^+l) < ^-|Je~^"^^e-^""'b-0 >* (B3) 



It is straightforward to show that 



. (p) • Je^'-'-'^e-'-^y'^ = l^e~'-'''^e-'''y'^J± (B4) 



with J± — J^; ± ijy the raising and lowering operators. Therefore 



ei(p) .X,„(p) = J^l±lD'*^,i^,e,0) (B5) 



and we have 



and 



Y«(p) = X,„(p) = ^l?2±l[DU^,9,O)e+{p)+Di^_,{^,0,O)e^{p)] (B6) 



Yg(p) = p X X,„(p) = -^^^^i±l[DCM()^0)e+{p) - DCi^p, 6, 0)e^{p)] (B7) 

Both y'j]1{p) and y'^-^\{p) are regular at poles — and 9 — n, so are Di^-^{if, 0, 0)e+(p) and D-l^_^{^, 0, 0)e_(p). 

In order to obtain the addition formula for Wigner D-function, we need to know the Euler angles of the product 
of two rotations with Euler angles {ip,9,0) and {ip\0',O). The composition rule of the transformation parameters of 
a Lie group is independent of the representations. For the rotation group, the J ~ 1/2 representation is the most 
convenient one to extract the composition rule. For J — 1/2, we have 

I?i(a,/3,7) = e-^'"^'"e-5'"^''e-^'^^T (B8) 
with (72,3 the Pauli matrices. We denote the combined rotation 

((p, 0, 0)D^- (^', 6*', 0) = L*^ (a, 7) (B9) 
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with a, (3, 7 the Euler angles of the combined rotation. Using Eq JBSp . we obtain 



cosfe-*("+^) 
sin|e3("-^) 



-sinfe-*("-r) 
cosfe*("+'^) 



cos I cos ^ei'^f'f"^ + sin | sin ^e~i(f-f'^ - cos f sin ^e^'^f'f'^ + sin | cos ^e'i 
-sin|cos^et(v-v') +cosf sin|^e-i(v'-v'') sin f sin |^ei(^-'^') 



■ cos ^ cos 2 



iv-v') 



The equivalence of each block on both sides gives rise to 

cos P = cos 9 cos 9' + sin 9 sin 6' cos{(p — (f') — p p' 

and 

g»(a+7) cos2^ =e+(p) •e_(p') 
ei("-7) sin2^ = -e+(p).e+(p') 



2 



-e_(p).e_(p') 



-("+7)cos2f =.e_(p).e+(p') 



(BIO) 

(Bll) 

(B12a) 
(B12b) 
(B12c) 
(B12d) 



These equations fix on the Euler angles of the combined rotation. Therefore, the composition rule of two rotations 
for representation of arbitrary J reads 



dfi(/3) cos2| 
e+(p) -e-lp') 



(B13) 
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